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^ \ ABSTRACT 

J> ■ We determine the abundance of i'-band drop-outs in the recently-released HSTj ACS 

£f~) ' Hubble Ultra Deep Field (UDF). Since the majority of these sources are likely to 

£N| , be z « 6 galaxies whose flux decrement between the F775W i'-band and F850LP z'- 

■ band arises from Lyman-alpha absorption, the number of detected candidates provides 
CO a valuable upper limit to the unextincted star formation rate at this redshift. We 

, demonstrate that the increased depth of UDF enables us to reach an 8 a limiting 

magnitude of z' AB = 28.5 (equivalent to 1.5 hj^ M yr _1 at z = 6.1, or 0.1 L* uv for 
the z w 3 [/-drop population), permitting us to address earlier ambiguities arising 
• from the unobserved form of the luminosity function. We identify 54 galaxies (and 

Qh only one star) at z' AB < 28.5 with (i! — z')ab > 1-3 over the deepest llarcmin 2 

q , portion of the UDF field. The characteristic luminosity (L*) is consistent with values 

■ observed at z « 3. The faint end slope (a) is less well constrained, but is consistent 
| with only modest evolution. The main change appears to be in the number density 

C$ . (*£*)• Specifically, and regardless of possible contamination from cool stars and lower 

k* ' redshift sources, the UDF data support our previous result that the star formation 

. £^ | rate at z « 6 was approximately x6 less than at z « 3 (Stanway, Bunker & McMahon 

^ ■ 2003). This declining comoving star formation rate (0.005 hm M yr 1 Mpc 3 at z w 6 

5-j ' at Ljjv > 0.1 L* for a salpeter IMF) poses an interesting challenge for models which 

suggest that Ljjv > 0.1 L* star forming galaxies at z ~ 6 reionized the universe. The 
short-fall in ionizing photons might be alleviated by galaxies fainter than our limit, or 
a radically different IMF. Alternatively, the bulk of reionization might have occurred 
at z ^> 6. 
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individual: SBM03#1 - galaxies: high redshift - ultraviolet: galaxies 



1 INTRODUCTION 

There has been considerable progress over the past decade in 
locating galaxies and QSOs at high redshifts. These sources 
have enabled us to probe the Universe at early epochs where 
its physical characteristics are fundamentally different from 
those at the present epoch. Observations of the most dis- 
tant z > 6.2 QSOs (Becker et al. 2001, Fan et al. 2002) 
show near-complete absorption at wavelengths shortward of 
Lyman-Q (Gunn & Peterson 1965), suggesting an optical 
depth in this line that implies a smooth neutral hydrogen 
fraction which is increasing rapidly with redshift at this 
epoch. Temperature-polarization cross-correlations in the 
cosmic microwave background from WMAP indicate that 



the Universe was completely neutral at redshifts of z > 10 
(Kogut et al. 2003). 

Although there is a growing consensus that cosmic 
reionization occurred in the redshift interval 6 < z < 15, a 
second key question is the nature of the sources responsible 
for this landmark event. Optical and X-ray studies to z ~ 6 
suggest the abundance of active galactic nuclei (AGN) at 
early epochs is insufficient when account is taken of the rel- 
evant unresolved backgrounds (Barger et al. 2003). A more 
promising source is star-forming galaxies whose early an- 
cestors may be small and numerous. Along with the escape 
fraction for the ionizing photons from the massive and short- 
lived OB stars in such sources, a major observational quest 
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in this respect is a determination of the global star formation 
rate at early epochs. 

In previous papers, our group has extended the Lyman- 
break technique (Steidel, Pettini & Hamilton 1995; Steidel 
et al. 1996) to address this question. Using the Advanced 
Camera for Surveys (ACS, Ford et al. 2002) on the Hubble 
Space Telescope (HST) with the sharp-sided SDSS F775W 
(?') and F850LP (z') filters, we located "i-drop" candidates 
with z' AB <25.6 at z ~ 6 for further study. In a series of 
papers, we have shown that this selection technique can ef- 
fectively locate z > 5.7 galaxies using ACS images from the 
HST Treasury "Great Observatory Origins Deep Survey" 
(GOODS; Giavalisco & Dickinson 2002). On the basis of 
GOODS-South photometric catalogues published by Stan- 
way, Bunker & McMahon (2003, hereafter Paper I), spectro- 
scopic follow-up using Keck/DEIMOS and Gemini/GMOS 
field demonstrated our ability to find high redshift galax- 
ies (Bunker et al. 2003, hereafter Paper II; Stanway et al. 
2004a). To address potential cosmic variance issues, we per- 
formed a similar analysis in the GOODS-North field, which 
yielded a consistent estimate of the surface density of z ~ 6 
star forming sources (Stanway et al. 2004b, hereafter Paper 
III). 

Although our initial study (Papers I-III) revealed the 
importance of ascertaining the difficult spectroscopic verifi- 
cations, and highlighted the problem of contamination from 
Galactic stars, we nonetheless determined that the abun- 
dance of confirmed star forming galaxies at z ~ 6 must be 
less than that expected on the basis of no evolution from the 
well-studied z ~ 3 — 4 Lyman break population (Steidel et 
al. 1999) . Working at the robustly-detected bright end of the 
luminosity function, in Paper I we showed that the comoving 
star formation density in galaxies with z' AB < 25.6 is « 6x 
less at z ~ 6 than at z ~ 3. Our z' AB < 25.6 flux limit corre- 
sponds to > 15 h^Q Mq yr _1 at z = 5.9, equivalent to Ly V 
at z rs 3. In Papers I-III we restricted our analysis to lumi- 
nous galaxies (where we take "luminous" to mean L > L* 
for the rest-UV). Other groups have claimed less dramatic 
evolution or even no evolution in the volume-averaged star 
formation rate, based on the same fields (Giavalisco et al. 
2004; Dickinson et al. 2004) and similar HST/ ACS data sets 
(Bouwens et al. 2003; Yan, Windhorst & Cohen 2003), but 
these groups work closer to the detection limit of the im- 
ages and introduce large completeness corrections for the 
faint source counts. The major uncertainty in converting 
the abundance of our spectroscopically-confirmed sample in 
the GOODS fields into a z ~ 6 comoving star formation rate 
is the form of the luminosity function for faint, unobserved, 
sources. As discussed in Paper III, if the faint end of the 
luminosity function at z ~ 6 was steeper than that at lower 
redshift, or if L* was significantly fainter, a non-evolving 
star formation history could perhaps still be retrieved. 

The public availability of the Hubble Ultra Deep Field 
(UDF; Beckwith, Somerville & Stiavelli 2003) enables us to 
address this outstanding uncertainty. By pushing the counts 
and the inferred luminosity function of i'-band drop-outs at 
2 w 6 to a limiting lower luminosity equivalent to one well 
below I/3 for the z w 3 population, it is possible to refine the 
integrated star formation rate at z ~ 6. In this paper we set 
out to undertake the first photometric analysis of i'-drops in 
the UDF. Our primary goal is to understand the abundance 
of fainter objects with characteristics equivalent to those of 



z ~ 6 sources and address uncertainties in the global star 
formation rate at this redshift. 

The structure of the paper is as follows. In Section 
we describe the imaging data, the construction of our cat- 
alogues and our i'-drop selection. In Section [3] we discuss 
the luminosity function of star-forming sources, likely con- 
tamination on the basis of earlier spectroscopic work, and 
estimate the density of star formation at z « 6. Our conclu- 
sions are presented in Section |3] Throughout we adopt the 
standard "concordance" cosmology of Qm — 0.3, S1a = 0.7, 
and use /170 = //o/70kms~ 1 Mpc" 1 . All magnitudes are on 
the AB system (Oke & Gunn 1983). 



2 HST IMAGING: OBSERVATIONS AND 
/-DROP SELECTION 

2.1 HST/ ACS Observations 

The Hubble Ultra Deep Field (UDF) is a public HST survey 
made possible by Cycle 12 STScI Director's Discretionary 
Time programme GO/DD-9978 executed over September 
2003 - January 2004. For the present program, the HST 
has imaged a single ACS Wide Field Camera (WFC) tile 
(11.5 arcmin 2 ) for 400 orbits in 4 broad-band filters (F435W 
B-band for 56 orbits; F606W U-band for 56 orbits; F775W 
i'-band for 144 orbits; F850LP z'-band for 144 orbits). The 
UDF field lies within the Chandra Deep Field South (CDF- 
S) with coordinates RA=03 fe 32 m 39!0, Decl --27°47'29'.'l 
(J2000). As the UDF represents the deepest set of images 
yet taken, significantly deeper than the /-band exposures of 
the Hubble Deep Fields (Williams et al. 1996; 1998), and 
adds the longer-wavelength z'-band, it is uniquely suited to 
the goals of our program. 

The WFC on ACS has a field of 202" x 202", and a 
pixel scale of 0.05". The UDF lies within the survey area of 
GOODS-South area (Giavalisco et al. 2004), surveyed using 
ACS with the same filters to shallower depth (3,2.5,2.5 & 5 
orbits in the B, V, i' & z' bands). The UDF was observed 
at two main orientations differing by 90 degrees, and within 
each of these data was taken in 2 blocks rotated by 4deg 
(orientations of 310,314,40 & 44 deg). A 4-point dither box 
spanning 0.3 arcsec was used, with half-pixel centres to im- 
prove the sampling. During each "visit" , there were 3 larger 

3 arcsec dithers to span the WFC inter-chip gap. 

For our analysis we use the reduced UDF data vl.O 
made public by STScI on 09 March 2004. The pipeline re- 
duction involved bias/dark current subtraction, flat-fielding, 
and the combination of background-subtracted frames re- 
jecting cosmic ray strikes and chip defects. The resulting re- 
duced images had been "drizzled" (Fruchter & Hook 2002) 
using the "MultiDrizzle" software (Koekemoer et al. 2004) 
on to a finer pixel scale of 0.03", to correct for geomet- 
ric distortion and to improve the sampling of the point 
spread function (PSF). The UDF data has been placed on 
the same astrometric system as the GOODSvl.O images of 
the UDF 1 . The photometric zeropoints adopted were those 
provided by STScI for the UDF vl.O data release: 25.673, 
26.486, 25.654 & 24.862 for the B, V, i' & z filters, where 
mag AS = zeropoint — 2.5 log 10 (counts/s). We have corrected 

1 Available from ftp://archive.stsci.edu/pub/hlsp/goods/vl 
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for the small amount of foreground Galactic extinction to- 
ward the CDFS using the CO&E/DIRBE & IRAS/ISSA 
dust maps of Schlegel, Finkbeiner & Davis (1998). The opti- 
cal reddening is E(B — V) = 0.008, equivalent to extinctions 24 - 
of A F775 = 0.017 & Afssolp = 0.012. 



2.2 Construction of Catalogues 

Candidate selection for all objects in the field was performed 
using version 2.3.2 of the SExtractor photometry package 
(Bertin & Arnouts 1996). As we are searching specifically 
for objects which are only securely detected in z', with min- 
imal flux in the i'-band, fixed circular apertures 0'.'5 in di- 
ameter were trained in the z'-image and the identified aper- 
tures used to measure the flux at the same spatial location 
in the i'-band image by running SExtractor in two-image 
mode. For object identification, we adopted a limit of at 
least 5 contiguous pixels above a threshold of 2a per pixel 
(0.0005 counts/pixel/s) on the data drizzled to a scale of 
/ /03 pixel -1 . This cut enabled us to detect all significant 
sources and a number of spurious detections close to the 
noise limit. As high redshift galaxies in the rest-UV are 
known to be compact (e.g., Ferguson et al. 2004, Bremer 
et al. 2004, Bouwens et al. 2004), we corrected the aper- 
ture magnitudes to approximate total magnitudes through 
a fixed aperture correction, determined from bright compact 
sources: —0.11 mag in i'-band and —0.14 mag in z'-band, the 
larger latter correction arising from the more extended PSF 
wings of the z'-band. 

The measured noise in the drizzled images underesti- 
mates the true noise as adjacent pixels are correlated. To 
assess the true detection limit and noise properties, we ex- 
amined the raw ACS/WFC images from the HST archive 
and measured the noise in statistically-independent pixels. 
For the 144-orbit z'-band, we determine that the 8 a detec- 
tion limit is z' AB = 28.5 for our O'.'5-diameter aperture. This 
is consistent with the noise decreasing as y/time from the 5- 
orbit GOODSvl.O to the 144-orbit UDF z'-band. We adopt 
this high S/N — 8 cut as our conservative sample limit. We 
trimmed the outermost edges where fewer frames overlapped 
in order to exploit the deepest UDF region, corresponding to 
a survey area of 11 arcmin 2 . From the output of SExtractor 
we created a sub-catalogue of all real objects brighter than 
z' AB < 28.5 mag (8a in a O'.'5-diameter aperture), of which 
63 appear to be promising i'-band dropouts (see §2.3) with 
(i' - z') AB > 1.3. 

To quantify possible incompleteness in this catalogue, 
we adopted two approaches. First we examined the recov- 
ery rate of artificial galaxies created with a range of total 
magnitudes and sizes. We used de Vaucouleurs r 1 ' 4 and ex- 
ponential disk profiles, convolved with the ACS/WFC PSF 
derived from unsaturated stars in the UDF images. Secondly 
we created fainter realisations of the brightest i'-dropout in 
the UDF confirmed to be at high redshift (SBM03#1 with 
z' AB = 25.4, confirmed spectroscopically to be at z = 5.83 
by Stanway et al. 2004b; Dickinson et al. 2004). By excising 
a small region around this i'-dropout, scaling the sub-image 
to a fainter magnitude, and adding it back into the UDF 
data at random locations, we assessed the recoverability as 
a function of brightness. For such objects we recover 98% 
of the simulated sources to z' AB = 28.5, the remainder be- 
ing mainly lost via source confusion through overlapping 
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Figure 1. The completeness (normalized to unity) for artificial 
galaxies added to the UDF z'-band image, as a function of total 
magnitude and half-light radius; we re-ran SExtractor on this 
image to assess the fraction of artificial galaxies recovered. The 
completeness is > 97% for R h < 0"2 and z' AB < 28.5. 

objects. From these analyses, we determine that, for un- 
resolved sources (r^ = 0'.'05), we are complete at our 8 a 
limit of z' AB = 28.5, and are 97% complete at this magni- 
tude for Vh = 0'.'2 (Figure 0. For objects with larger half- 
light radii we will underestimate the z'-band flux due to our 
O'.'5-diameter photometric aperture. However, this effect is 
small for our sample of compact sources (Table 1 lists both 
the O'.'5-diameter magnitudes with an aperture correction 
which we adopt, and the SExtractor "MAG_AUTO" esti- 
mate of the total magnitude using a curve-of-growth: these 
are broadly consistent). 

At the relatively bright cut of z' AB < 25.6 used in Pa- 
per I from the GOODSv0.5 individual epochs, the UDF data 
is 98% complete for sources as extended as ru = 0.5 arcsec. 
Interestingly, we detect no extended (low surface brightness) 
i'-drops to this magnitude limit in addition to SBM03#1 
(Papers I, III) in the deeper UDF data. This supports our 
assertion (Paper I) that the i'-drop population is predom- 
inantly compact and there cannot be a large completeness 
correction arising from extended objects (c.f. Lanzetta et al. 
2002). The ACS imaging is of course picking out Hn star 
forming regions, and these UV-bright knots of star forma- 
tion are typically < 1 kpc (< 0'.'2 at z ~ 6) even within large 
galaxies at low redshift. 

2.3 z ~ 6 Candidate Selection 

In order to select z ~ 6 galaxies, we use the Lyman break 
technique pioneered at z ~ 3 using ground-based telescopes 
by Steidel and co-workers and using HST by Madau et al. 
(1996). At z ~ 3 — 4 the technique involves the use of three 
filters: one below the Lyman limit (A rC st = 912 A), one in the 
Lyman forest region and a third longward of the Lyman- 
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a line (A rcst = 1216 A). At z « 6, we can efficiently use 
only two niters, since the integrated optical depth of the 
Lyman-Q forest is 3> 1 (see Figure|5J rendering the shortest- 
wavelength filter below the Lyman limit redundant. The key 
issue is to work at a sufficiently-high signal-to-noise ratio 
that i'-band drop-outs can be safely identified through de- 
tection in a single redder band (i.e., SDSS-2'). This approach 
has been demonstrated to be effective by the SDSS collabo- 
ration in the detection of z ~ 6 quasars using the i'- and 
z'-bands alone (Fan et al. 2001). The sharp sides of the 
SDSS filters assist in the clean selection using the photo- 
metric redshift technique. In Figures |21&Bl we illustrate how 
a colour cut of (j! — z')ab > 1-5 (used in Papers I-III) can 
be effective in selecting sources with z > 5.7. Here we re- 
lax this cut to (i' — z')ab > 1.3 to recover most galaxies at 
redshifts z > 5.6, but at the expense of potentially larger 
contamination by z ~ 1 — 2 ellipticals. Near-infrared colours 
from the NICMOS imaging of the UDF should identify these 
Extremely Red Objects (EROs), and we consider this in a 
companion paper (Stanway, McMahon & Bunker 2004c). 
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Six of the 63 candidate i'-dropouts in our z' AB < 28.5 
UDF catalogue were identified visually as different regions 
of the same extended source, and where these were within 
our aperture diameter of 0'.'5 the duplicates were eliminated 
from the final selection. One spurious i'-drop arose from the 
diffraction spikes of bright stars due to the more extended 
PSF in the z'-band compared with that in the i'-band. Only 
one of the i'-dropouts is unresolved (Figure [SJ. This is the 
brightest at z' AB — 25.3 (#11337 in Table 1), detected in the 
U-band image and removed from our catalogue of potential 
z « 6 objects as a probable star. At the edge of the UDF 
frame (and outside the central llarcmin 2 region of lowest 
noise where we do our main analysis) there is a second un- 
resolved i'-drop with z' AB = 25.2, first identified in Paper I 
(SBM03#5), where we argued that the near-IR colours are 
likely to be stellar. It is interesting that the level of stellar 
contamination in the UDF i'-drops is only 2%, compared 
with about one in three at the bright end (z' AB < 25.6, Pa- 
pers I & III). This may be because we are seeing through the 
Galactic disk at these faint magnitudes to a regime where 
there are no stars at these faint limiting magnitudes. 



From our original list of 63 i'-drops, 6 duplications 
were removed, along with one diffraction spike artifact. 
The remaining objects satisfying our (i' — z') AB > 1.3 & 
z' AB < 28.5 selection criteria are detailed in Table 1, of which 
54 are good candidate z « 6 galaxies, along with the proba- 
ble star #11337, and another objected (#46574) detected in 
U-band. The surface density of i'-drops as a function of lim- 
iting magnitude is shown in Figure |7| None of the i'-drops 
(with the exception of the Galactic star) are detected in the 
B-band image of the UDF, to a 3 a limit of B AB > 29.2 in a 
O'.'5-diameter aperture, as would be expected for the z ~ 6 
interpretation where the B-band filter covers wavelengths 
below the 912 A Lyman limit. 



Figure 2. The ACS-i' and -z' bandpasses overplotted on the 
spectrum of a generic z = 6 galaxy (solid line), illustrating the 
utility of our two-filter technique for locating 2^6 sources. 



3 SELECTION EFFECTS AND THE 
LUMINOSITY FUNCTION OF STAR 
FORMING GALAXIES AT2~6 

3.1 Estimate of Star Formation Rate from the 
Rest-UV 

We will base our measurement of the star formation rate for 
each candidate on the rest-frame UV continuum, redshifted 
into the z'-band at z ~ 6 and measured from the counts in a 
O'.'5-diameter aperture (with an aperture correction to total 
magnitudes, Section ^. 2H . In the absence of dust obscuration, 
the relation between the flux density in the rest-UV around 
« 1500 A and the star formation rate (SFR in MQyr -1 ) 
is given by Luv = 8 x 10 2r SFRergs s^ 1 Hz~* from Madau, 
Pozzetti & Dickinson (1998) for a Salpeter (1955) stellar 
initial mass function (IMF) with 0.1 M < M* < 125 M . 
This is comparable to the relation derived from the mod- 
els of Leitherer & Heckman (1995) and Kennicutt (1998). 
However, if a Scalo (1986) IMF is used, the inferred star 
formation rates will be a factor of w 2.5 higher for a similar 
mass range. 

Recognising the limitations of our earlier studies (Pa- 
pers I-III) which by necessity focussed on the brighter i'- 
drops, we now attempt to recover the z ~ 6 rest-frame UV 
luminosity function from the observed number counts of i'- 
drops to faint magnitudes in the UDF. Although our colour 
cut selects galaxies with redshifts in the range 5.6 < z < 7.0, 
an increasing fraction of the z'-band flux is attenuated by 
the redshifted Lyman-a forest. At higher redshifts we probe 
increasingly shortward of A rcst = 1500 A (where the lumi- 
nosity function is calculated) so the fe-corrections become 
significant beyond z m 6.5. 

Figure |5] demonstrates this bias and shows the limit- 
ing star formation rate as a function of redshift calculated 
by accounting for the filter transmissions and the blanket- 
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Table 1. i'-band dropouts in the UDF. The two stars are above the line - all others are spatially resolved. Our ID and the corresponding 
match from the UDF catalogues released by STScI are listed. Where two close i'-drops lie within our 0"5-diamcter aperture, the flux 
only counted once in the star formation total — those IDs and star formation rates in parentheses are not counted. The star formation 
rates assume the i'-drops lie at z = 6.0, the expected median redshift of our sample. The z' AB (total) is the SExtractor "MAG-AUTO". 
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5 double. * star SBM03# 5 (Paper I), outside central UDF. 1 SBM03#1 (Paper I); SiD002 (Dickinson et al. 2004). 2 SiD025 (Dickinson 
et al. 2004). 3 46574 has a close neighbour visible in the u-band (i.e. low redshift.) 4 46234 is close to 46223. 5 27029 is close to 27032. 
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Figure 3. Model colour-redshift tracks for galaxies with non- 
evolving stellar populations (from Coleman, Wu & Weedman 1980 
template spectra). The contaminating 'hump' in the («' — z') 
colour at 2 rj 1 - 2 arises when the Balmer break and/or the 
4000 A break redshifts beyond the i'-filter. 



ting effect of the intervening Lyman-a forest. By intro- 
ducing the small fe-correction to A rcst = 1500 A from the 
observed rest-wavelengths longward of Lyman-a redshifted 
into the z'-band we can correct for this effect. We consid- 
ered a spectral slope of (3 = —2.0 (where fx oc A' 3 ) appro- 
priate for an unobscured starburst (flat in /„), and also a 
redder slope of f3 — —1.1 which appropriate for mean red- 
dening of the z ~ 3 Lyman break galaxies given by Meurer 
et al. (1997). A more recent determination for this popu- 
lation by Adelberger & Steidel (2000) gives /9 = —1.5, in 
the middle of the range. At our 8<r limiting magnitude of 
z 'ab = 28.5, we deduce we can detect unobscured star for- 
mation rates as low as 1.0 [1.1] Mq yr _1 at 5.6 < z < 5.8 
and 1.5 [1.7] /if 2 Mq yr _1 at z < 6.1 for spectral slope 
P= -2.0 [-1.1] (Figure^. 

Recognising that contamination by interlopers will only 
reduce the value, we now compare the comoving star for- 
mation rate deduced for z ~ 6 galaxies based on our can- 
didate i'-dropout source counts with predictions based on 
a range of rest-frame UV luminosity functions. For conve- 
nience we assume that there is no evolution over the sam- 
pled redshift range, 5.6 < z < 6.5, spanned by the UDF 
data (equivalent to a range between 0.8 — 1.0 hj^ Gyr after 
the Big Bang). We take as a starting point the luminos- 
ity function for the well-studied Lyman-break [/-dropout 
population, reported in Steidel et al. (1999), which has a 
characteristic rest-UV luminosity m* R — 24.48 (equivalent 
to M£(150QA) = -21.1 mag or L* 3 = 15 M© yr" 1 for 
our cosmology). The faint end slope of the Schechter func- 
tion at z ~ 3 is relatively steep (a = —1.6) compared with 
a — —1.0 to —1.3 for lower-redshift galaxy samples (e.g., 
Lilly et al. 1995; Efstathiou et al. 1988; Blanton et al. 2003 - 
see Gabasch et al. 2004 for recent determinations at 1500 A). 
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Figure 4. Colour-magnitude diagram for the UDF data with the 
limit z' AB < 28.5 and (i' — z')j\g = 1.3 colour cut shown (dashed 
lines). As discussed in the text, such a catalogue could be con- 
taminated by cool stars, EROs and wrongly identified extended 
objects and diffraction spikes but nonetheless provides a secure 
upper limit to the abundance of z?s 6 star forming galaxies. Cir- 
cles and arrows (lower limits) indicate our i'-drop candidate 2 fs 6 
galaxies. The solid circle is the spectroscopically-confirmed galaxy 
SBM03#1 (Stanway et al. 2004b; Dickinson ct al. 2004), and the 
asterisk is the only unresolved i'-drop in our UDF sample, the 
probable star #11337. 



The characteristic comoving number density at z 

l3 AJ„„-3„ -1 



3 is 



We adopt two approaches to determining the galaxy 
number density and star formation density at z w 6: the first 
method (Section l3.1.1H is the one used in Papers I & III, an 
application of the "effective volume" technique (Steidel et al. 
1999). The second method (Section l3.1.3t involved compar- 
ing the measured surface density of i'-dropout z ~ 6 galaxies 
with that predicted on the assumption they have the same 
characteristics as the [/-dropout population at z ~ 3. 



3.1.1 Effective Survey Volume 

We have followed the approach of Steidel et al. in calculating 
the effect of luminosity bias on our sample of z ~ 6 LBGs. 
We account for the fc-correction: as redshift increases, the 
z'-band samples light in the rest-frame of the galaxies at 
wavelengths that are increasing ly far to the blue of 1500 A, 
where the LBGs' luminosity function was calculated. Addi- 
tionally, at redshifts z > 6, Lyman-a absorption from the 
forest enters the z'-band and makes galaxies fainter still, as 
there is incomplete coverage of the filter by the continuum 
longward of Lyman-a. Accounting for these luminosity and 
redshift biases, we compute an effective survey volume using 



Ves (m) 
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Figure 5. The distribution of angular sizes (half-light radius, 
Rh, in arcseconds) for objects in our z'-band selected catalogue. 
Our i'-drop candidate z PS 6 are marked as open circles, with the 
confirmed z = 5.8 galaxy SBM03#1 a solid circle. The i'-drops 
appear to be compact but resolved (the stellar locus at 0'.'05 is 
clearly visible). The asterisk denotes the only unresolved i'-drop 
in our UDF sample, the probable star #11337. 



where p(m, z) is the probability of detecting a galaxy at 
redshift z and apparent z magnitude m, and dz 4*- is the 
comoving volume per unit solid angle in a slice dz at redshift 
z. We integrate over the magnitude range we are sensitive 
to, and over the redshift range 5.6 < z < 7.0 from our colour 
selection, and calculate that for a spectral slope of f3 = —2.0 
(i.e., flat in /„) the effective comoving volume is 40 per cent 
the total volume in the range 5.6 < z < 7.0 (i.e., the same 
as 5.6 < z < 6.1). For our 11 arcmin 2 survey area (excluding 
the edge regions of the UDF where fewer frames overlap) 
this is a comoving volume of 2.6 x 10 4 Mpc 3 . 

Hence we calculate a volume-averaged comoving 
star formation density at z « 6 of (0.005 ± 
0.001) /i7oM yr _1 Mpc -3 for the w 50 i'-dropout galax- 
ies with z'(AB) < 28.5 (L uv > 0.1 L%). This is plotted 
on the Madau-Lilly diagram (Figure UTt . Data from other 
groups are shown on this figure, where we have corrected 
all the datasets to the same limiting star formation rate 
of 1.5 h^Q Mq yr _1 (i.e., typically integrating their claimed 
luminosity functions down to « 0.1 L*s) to provide a fair 
comparison of evolution. Integrating the luminosity func- 
tion down to ~ 0.1 L*, as here, represents most of the to- 
tal luminosity density for faint end slopes a > —1.6 (com- 
pared with integrating to zero luminosity). If we assume 
that the Schechter function holds for the unobserved faint 
galaxies with L < 0.1 L*3, then the observed population 
with L > 0.1 Lg represents (87.5%, 78.9%, 56.4%, 32.4%, 
17.4%) of the total star formation rate for faint-end slopes 
a = (-1.1, -1.3, -1.6, -1.8, -1.9). 
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Figure 6. Limiting star formation rate as a function of redshift 



for the UDF catalogue with 



< 28.5 mag (8<r). Star forma- 



tion rates are inferred from the rest-frame 1500 A flux (Madau, 
Pozzetti & Dickinson 1998) taking account of fc-corrections, filter 
transmission and blanketting by Lyman-a absorption. The solid 
line assumes a spectral slope /3 = —2.0 (where fx oc A* 3 ) appropri- 
ate for an unobscured star burst, and the dotted line has = —1.1 
(corresponding to mean reddening of z «3 3 Lyman break galax- 
ies given in Meurer et al. 1997). The limit as a fraction of L* A 



(L*[1500A] at z « 3, equivalent to SFR{ 



15/i 7£ ;M Q yr- 



from Steidel et al. 1999) is shown on the right axis. Our colour 
selection should remove most z < 5.6 galaxies (solid vertical line), 
and our average i'-drop redshift for z' < 28.5 should be z m 6.0 
(vertical dot-dash line): we are sensitive as faint at 0.1 Lg at this 
redshift. 



3.1.2 Colour Selection: spectral slope and forest 
blanketting 

We model the effect of the break below the Lyman-a emis- 
sion line due to blanketting by the forest, where the contin- 
uum break Da (Oke & Korycansky 1982) is defined as 



D A = 



/„(1050 - 1170 A) c 
/„(1050- 1170 A) 



prcd 



(i) 



We assumed flux decrements of Da = 0.9 — 1.0, consistent 
with that observed in the z > 5.8 SDSS QSOs (Fan et al. 
2001). We find that lowering Da reduces the completeness 
in the lowest redshift bin 5.6 < z < 5.8 for a (i' — z')ab > 1-5 
colour cut. A (i' — z')ab > 1.3 cut improves the selection 
somewhat but at the risk of higher contamination from red 
objects at z ~ 1 — 2: we consider this in Stanway, McMahon 
& Bunker (2004c). 

We find that altering the spectral slope /3 of the i'- 
drop spectral energy distribution (fx tx A - * 3 ) over the range 
— 1.1 > fi > —2.0 changes the predicted number of i! - 
dropouts by only « 10%. 
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3.1.3 Surface Density Predictions 

First, we compare our observed number of i'-band dropout 
galaxies with a simple no-evolution model, which assumes 
the same luminosity function for Lyman break galaxies at 
z — 6 as at z = 3 (with a faint-end slope a = —1.6, spectral 
slope [3 — —2.0 and Lyman-a forest decrement D A = 1.0). 
This no-evolution model would predict 169 galaxies satisfy- 
ing our i' AB < 28.5 & (i' — z')ab > 1.3 selection with a total 
star formation rate of 866 /if 2 Mq yr -1 . This compares with 
our observed number of 54 i'-drops (l/3rd that predicted), 
which have a total star formation rate of 140 /if 2 Mq yr -1 
(l/6th of the no-evolution prediction). The predicted me- 
dian redshift of our i'-drop sample for the no-evolution 
model is z = 5.95, with the luminosity- weighted average 
2 = 6.05. 

Clearly, evolution in the UV luminosity function of Ly- 
man break galaxies is required. To fit this, we constructed 
a grid of models based upon the z w 3 luminosity function, 
varying a between —1.1 and —1.9, and L* between 0.3 L3 
and 2 L3. We leave the normalization of the luminosity func- 
tion, <&*, as a free parameter, 

We minimize % 2 f° r °ur grid of model luminosity func- 
tions: our best fit (Figure[HJ is compatible with no evolution 
of L* from 2 w 3, but a large decline in the comoving space 
density, <!>* (by about a factor 6 relative to 2 ~ 3). The faint 
end slope is less well constrained, although no evolution is 
compatible with the results. At the faintest magnitude bin, 
there modestly higher counts, perhaps indicating a slightly 
steeper a if the results at the faintest magnitudes are to be 
trusted (Figure |5J. 

3.1.4 Comparison with i -drop number counts from other 
groups 

Recognizing the very limited area of the UDF and the prob- 
lems of cosmic variance, it is nonetheless interesting to com- 
pare our measured i'-drop number counts with previous 
determinations from shallower data sets. The surface den- 
sity derived in Paper I to z' AB — 25.6 is consistent with 
the present data - we detect only one resolved i'-dropout 
this bright: SBM03#1. Note that the UDF pointing was se- 
lected to include this object. No other spatially-resolved i'- 
dropouts are detected to z' AB < 26.5, implying a surface den- 
sity of 0.1±0.1 arcmin -2 . This is in contrast with the density 
of 0.4 arcmin -2 measured by Giavalisco et al. (2004) to the 
same z' AB < 26.5 limit, and 0.5 ± 0.2 arcmin -2 from the 
completeness-corrected estimate of Bouwens et al. (2003) 2 , 

2 Note added in proof: a recent paper by Bouwens et al. (2004), 
based on number counts of i'-drops in the ACS parallel observa- 
tions to the NICMOS UDF field, significantly revises their previ- 
ous estimate of the number density of z' AB < 26.5 «'-drops from 
0.5 ±0.2 to 0.2 ±0.1 arcmin -2 (4 objects in 21 arcmin 2 , consistent 
with our UDF work presented here). The conclusion of Bouwens 
ct al. (2003) —that the comoving star formation at z ft: 6 is con- 
sistent with no evolution from z ft; 4- is revised in Bouwens et al. 
(2004) to be a factor of 3 decline from z = 3.8 to z ta 6. Using the 
evolution in comoving number density of (1 + z) -2,8 suggested 
by Bouwens et al. (2004), this fall in star formation rate at z = 6 
is consistent with our result of a factor of 6 decline from z = 3.0 
to z Ri 6 from the GOODS data in Stanway, Bunker & McMahon 
(2003), confirmed in this paper from the deeper UDF data. 
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Figure 7. Cumulative source counts per arcmin 2 of i'-dropouts 
as a function of z'-band magnitude. The new UDF data (over 
a smaller area of 11 arcmin 2 for z' Ag > 27.0) is compared with 
z' AB < 25.6 single epoch GOODSv0.5 ACS/WFC imaging over 
300 arcmin 2 (Papers I-III) and combined 5 epoch GOODSvl. 
images to z' AB < 27.0 (Stanway 2004). 

and the even higher surface density of 2.3 arcmin -2 (after re- 
moving stellar contamination) claimed by Yan, Windhorst 
& Cohen (2003), after correcting for a factor of 4 error in 
their original flux calibration (see Yan & Windhorst 2004). 
Clearly, there are large discrepancies from the various groups 
in the number density measured to the same limiting mag- 
nitude of z' AB < 26.5, with measurements up to a factor of 
20 higher than our UDF measurement (Yan & Windhorst 
2004). These discrepancies may be due to cosmic variance, 
or too many spurious sources in the samples of these teams, 
due to working close to the sensitivity limits. By using a high 
signal-to-noise (S/N — 8) cut, we guard against the low-S/N 
bias: where there are many more objects with intrinsically 
bluer colours that scatter up into our (i' — z') AB > 1.3 se- 
lection than there are real i-drops which scatter out of the 
colour selection through photometric errors. 

From Somerville et al. (2004) we estimate that the cos- 
mic variance for the UDF is 40%, assuming the 2 = 6 LBGs 
are clustered in the same way as the 2 = 3 LBGs and as- 
suming a volume of derived by scaling our UDF area with 
our wider-area GOODS data (with an effective volume of 
1.8 x 10 5 ftfo 3 Mpc 3 for the 146 arcmin 2 of GOODS-S, Pa- 
per I). Indeed, the spatial distribution of our i'-drops on 
the sky does indicate some clustering (Figure [TOt . and we 
had already flagged 6 of our candidates as being "double" 
sources (Table 1), with another 2 having near neighbours. 
In the GLARE GMOS/Gemini spectroscopy of the GOODS- 
South i'-dropouts, Stanway et al. (2004a) have already spec- 
troscopically identified an overdensity at 2 = 5.8. 

3.2 Implications for Reionization 

The increased depth of the UDF enables us to resolve the 
uncertainties associated with the unobserved portion of the 
luminosity function (LF) for 2 w 6 sources. Our best-fit 
LF suggest little or no change in L* over 3 < 2 < 6, with 
a less well constrained but consistent with modest evolu- 
tion, implying the major evolution is a decline in space den- 
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Figure 8. Cumulative source counts per arcmin 2 of i'-dropouts 
as a function of 2 ; -band magnitude, with various normalisations 
of the characteristic number density at z xs 6, $g (in terms of the 
value at z rs 3, assuming Lg = and the same a as the 

z 3 Lyman break population (a = —1.6). Our faintest point 
:' AB = 26.5) is excluded from the fit 
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Figure 9. Cumulative source counts per arcmin 2 of i'-dropout as 
a function of z'-band magnitude, with various values of the faint 
end slope (a) assuming Lg = Lg and $g = $3/6. Our faintest 
point from the GOODS data (at z' AB = 26.5) is excluded from 
the fit due to incompleteness. 



sity (and global star formation rate) by ~ x6 at 2 Ri 6. 
This sharp decline, which must represent a lower limit to 
the true decline given the likelihood of contamination from 
foreground sources, suggests it may be difficult for luminous 
star-forming z ~ 6 i-dropout galaxies to be the main source 
of ionizing photons of the Universe. 

We attempt to quantify this by comparing with the es- 
timate of Madau, Haardt & Rees (1999) for the density of 
star formation required for reionization (their equation 27): 



PSFR : 



0.013 Meyr' 1 Mpc~ 

/esc 



0.08 



(I) <» 



This relation is based on the same Salpeter Initial Mass 
Function as we have used in deriving our volume-averaged 
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Figure 10. The spatial distribution of our UDF i'-drops on the 
sky (diamonds). The location of the confirmed z = 5.8 source from 
Paper I is marked (#1) as are two other sources just outside the 
UDF, spectroscopically identified at 2 = 5.8 — 5.9 by Stanway et 
al. (2004a). 



star formation rate. C is the concentration factor of neutral 
hydrogen, C = (phi) (pHi)~ 2 . Simulations suggest C w 30 
(Gnedin & Ostriker 1997). Our comoving star formation rate 
of 0.005 hro M Q yr _1 Mpc~ 3 from the i'-drop galaxies we de- 
tect is a factor of > 2.5 lower than the original Madau, 
Haardt & Rees (1999) requirement at z ~ 5. We have up- 
dated their equation 27 for the more recent concordance 
cosmology estimate of the baryon density of Spergel et al. 



(2003), Q b = 0.0224 



0.0457 h 70 , and for the predicted 



mean redshift of our sample (z = 6.0): 



PSFR : 



0.026 Meyr" 1 Mpc~ 

/esc 



(l±±Y ( »6 fero V (C_\ 
\ 7 J { 0.0457/ V30/ 



(3) 



The escape fraction of ionizing photons (/ eS c) for high- 
redshift galaxies is highly uncertain (e.g., Steidel, Pettini & 
Adelberger 2001), but even if we take / OS c = 1 (no absorption 
by H 1) this estimate of the star formation density required 
is a factor of w 5 higher than our measured star formation 
density of 0.005 hjo Mq yr _1 Mpc~ 3 at z w 6 from galax- 
ies in the UDF with SFRs > 1.5 /i TO 2 Mq yr~ x . For faint end 
slopes of a — 1.8 — » —1.3 galaxies with L > 0.1 L* account 
for 32 — 80% of the total luminosity, so would fall short of 
the required density of Lyman continuum photons required 
to reionize the Universe. If the faint-end slope is as steep as 
a w —1.9 then there would just be enough UV Lyman con- 
tinuum photons generated in star forming galaxies at 2 w 6 
(assuming a Salpeter IMF), but the requried escape fraction 
for complete reionization would still have to be implausibly 
high (/esc ~ 1, whereas all high-z measurements to date in- 
dicate /e S c -C 0.5: Fernanadez-Soto, Lanzetta & Chen 2003; 
Steidel, Adelberger & Pettini 2001). 

AGN are also under-abundant at these epochs (e.g., Di- 
jstra, Haiman & Loeb 2004). If star forming galaxies at red- 
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Figure 11. An updated version of the 'Madau-Lilly' diagram 
(Madau et al. 1996; Lilly et al. 1996) illustrating the evolution 
of the comoving volume-averaged star formation rate. Our work 
from the UDF data is plotted a solid symbol. Other determina- 
tions have been recalculated for our cosmology and limiting UV 



luminosity of 1.5 /i 70 Mq yr at z = 6.1 (equivalent to 0.1 Lg 
at z « 3 from Steidel et al. 1999), assuming a slope of a = —1.6 
for 2 > 2 and a = -1.3 for z < 2. Data from the CFRS sur- 
vey of Lilly et al. (1996) are shown as open circles; data from 
Connolly et al. (1997) are squares; and the Lyman break galaxy 
work of Steidel et al. (1999) is plotted as crosses, of Fontana et 
al. (2002) as inverted triangles and that by Iwata et al. (2003) as 
an open diamond. Pentagons are from Bouwens, Broadhurst & 
Illingworth (2003). The upright triangles are the GOODS i'-drop 
results from Giavalisco et al. (2004). The three ACS estimates of 
Bouwens et al. (2003) are shown by small crossed circles and in- 
dicate three different completeness corrections for one sample of 
objects - the larger symbol is the recent re-determination using 
a new catalogue by this group from a deeper dataset (the UDF 
flanking fields - Bouwens et al. 2004); we have recomputed the 
comoving number density from the Bouwens et al. (2004) because 
of a discrepancy on the scale of their plot of star formation history 
(their Fig. 4 in|astro-ph7u403167 vl fc v2). 



shifts close to z — 6 were responsible for the bulk of reion- 
ization, then a very different initial mass function would be 
required, or the calculations of the clumping factor of neu- 
tral gas would have to be significantly over-estimated. Al- 
ternatively another low- luminosity population (e.g., forming 
globular clusters; Ricotti 2002) could be invoked to pro- 
vide some of the shortfall in ionizing photons. It is also 
plausible that the bulk of reionization occured at redshifts 
well beyond z = 6: the WMAP polarization data indicate 
Zreion > 10 (Kogut et al. 2003), and it is possible that the 
Gunn-Peterson troughs seen in z > 6.2 QSOs (Becker et al. 
2001; Fan et al. 2002) mark the very last period of a neutral 
IGM. 



4 CONCLUSIONS 

We summarize our main conclusions as follows: 

(i) We present an i'-dropout catalogue of z ~ 6 star form- 
ing galaxy candidates in the Ultra Deep Field (UDF) to a 
limiting flux (8<r) of z' AB < 28.5. This represents a substan- 



tial advance over the depths achieved in the GOODS cata- 
logues and enables us, for the first time, to address questions 
concerning the contribution of the faint end of the luminos- 
ity function. 

(ii) We detect 54 resolved sources with (i' — z')ab > 1.3 
in the deepest llarcmin 2 portion of the UDF and consider 
this to be an upper limit to the abundance of star forming 
galaxies at z ~ 6. 

(iii) Using simulations based on lower redshift data, we 
deduce that, regardless of contamination by foreground in- 
terlopers, the abundance of i'-dropouts detected is signifi- 
cantly less than predicted on the basis of no evolution in 
the comoving star formation rate from z — 3 to z — 6 (in- 
tegrating to the same luminosity level). The UDF data sup- 
ports our previous suggestions that the star formation rate 
at z ~ 6 was about x 6 less than at z ~ 3 (Stanway, Bunker 
& McMahon 2003). 

(iv) The inferred comoving star formation rate of 
0.005 h 7 o M e yr -1 Mpc -3 from L > 0.1 L* vv galaxies at 
z ~ 6 may poses a significant challenge for models which 
require that luminous star forming galaxies in the redshift 
range 6< z <10 are responsible for reionizing the Universe. 

(v) The contamination of our i'-drop sample of candidate 
z ss 6 galaxies by cool Galactic stars appears to be minimal 
at z' AB > 26, possibly because we are seeing beyond the 
Galactic disk at the faint magnitudes probed by the UDF. 

Note Added in Proof 



A recent preprint by Yan & Windhorst (astro-ph/0407493 ) 
independently repeats our selection of candidate z ~ 6 
galaxies in the Hubble Ultra Deep Field with (i 1 — z')ab > 
1.3 ( astro-ph/0403223 and this paper). The Yan & Wind- 
horst catalogue also pushes to fainter magnitudes than our 
z' AB < 28.5 limit, where the completeness corrections be- 
come significant. This subsequent independent analysis re- 
covers almost all of our original i'-band drop-out galaxies, 
and the catalogues agree at the 98% level (one discrepant 
object out of 50). In astro-ph/0407562 (Bunker & Stan- 
way 2004) we present a matched catalogue of these i-band 
dropouts in the Hubble Ultra Deep Field. 
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